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Protein Structure and Function 
by Comparative Model Building 

JONATHAN GREER 

Physical Biochemistry Laboratory 
Computer-Assisted Molecular Design 
Abbott Laboratories 
Abbott Park, North Chicago. Illinois 60064 

INTRODUCTION 

Comparative modeling methods have been used to extend the experimen- 
tally determined three-dimensional structures of proteins to new molecules 
whose structure is closely related. Such techniques have been applied to de- 
rive model structures for a-lactalbumin,' a-lytic protease,^ Streptomyces 
trypsin-like protein,^ Ca** binding proteins/ haptoglobin/ serine proteases/ 
including blood clotting factor Xa^ and very recently renin/ a member of 
the acid protease family. Thus, comparative model building has been widely 
used for producing tentative structures of biologically interesting and im- 
portant molecules. 

In this study, we employ comparative modeling methods to begin ex- 
ploring the nature of specificity between enzymes and their particular sub- 
strates. Initial structures for several serine proteases are derived.*"^ Suitable 
peptides from the known macromolecular substrates of these enzymes are 
also modeled onto the enzyme active site^ and their properties in relation 
to their respective enzymes are analyzed and compared. The ultimate goal 
is to gain a deeper and more detailed understanding of the molecular basis 
of enzyme-substrate and protein-ligand recognition and specificity. 



METHODS 

Comparative Modeling Method 

A detailed description of the comparative modeling methods for the serine 
proteases used in this work has been published.* Briefly, by comparing the 
experimentally known structures of chymotrypsin,'''" trypsin,'*-^^ and 
elastase,'^ ''' the three-dimensional structures of the serine proteases can be 
parsed into structurally conserved regions (SCRs) and variable regions (VRs) 
(see Fig. 1). Sequence homology among these proteins lies almost exclusively 
in the SCRs (see Figure 2 of ref. 6). 

To model a "new" serine protease, the sequence is aligned using the strong 
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homology in the SCRs. The coordinates for the main chain in the SCRs are 
taken from any one of the known structures (since they are almost identical, 
see Figure 1 of ref. 6). The side chains are "mutated" to fit the sequence of 
the new protein. 

Modeling the VRs is much more challenging. For each VR, the various 
conformations found amongst the known structures (Fig. 1) are examined 
as to length (see Table 3 of ref. 6) and residue character. If one of the known 
conformations fits the new sequence, its main chain coordinates are used 
directly with suitable replacement of side chains. Otherwise, modehng using 
energetics is necessary to achieve a reasonable tentative conformation for 
the respective VR. Modeling studies on eight serine protease sequences'*-' show 
that more than 50% of the VRs can be modeled directly from the known 
structures whereas less than 5% fall in a class of large additions where 
modeling is not practical in the near future. 



Modeling the Serine Proteases and Their Substrates 

Table 1 lists the substrates and enzymes used in this study. The two serine 
proteases, blood clotting factors Xa and IXa, were sequenced by Titani et 
al.'' and Katayama et al.,'^ respectively. The sequence alignments used for 
these two proteins are given in Figure 3 of ref. 6. The SCRs were modeled 
from elastase and the VRs as shown in Table 2. In the case of factor Xa 
(see ref. 7) three loops cannot be built without further energy analysis; how- 
ever, these three VRs are distant from the active site region and thus do not 
concern us in this study. Factor IXa is closely related but not identical to 
factor Xa in the size of its VRs as can be seen from Figure 3 and Table 3 
of ref. 6 and from Table 2. In factor IXa, only two loops have no direct 
known structure as a model. One of them, the VR at 36-38, lies just on the 
border of the substrate binding pocket. The other is distant from the active 
site region, as it was in factor Xa. 

The substrate peptides, consisting of residues P4,..., Pi, Pi',..., P/ (using 
the standard substrate nomenclature"), were modeled after the crystal- 
lographically determined conformations of bovine pancreatic trypsin inhib- 
itor,'i of soybean trypsin inhibitor, and of di- and tripeptide chloromethyl 
ketones bound to y-chymotrypsin^"-^' and to Streptomyces griseus protease 
B." All these inhibitors have a common main chain conformation between 
residues P3 and P3' (see ref. 7 for a more complete discussion of this). The 
side chain positions for residues P3 and P3' as well as the position of residue 



TABLE 1. Substrates and Enzymes Used 



Substrate 


Enzyme 


prothrombin 


factor Xa 


factor X 


factor IXa 


trypsinogen 


trypsin 
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TABLE 2. Model Structures for the VRs in Factors Xa and IXa 



V K 


iFactor Xa 


Factor IXa 


Model 


Residues built 


Model 


Residues built 


23-25 


elastase 




elastase 





36-38 


rhvmntrvnsin 

^11 Jr Lit \J Lt y L/Olli 


33-40 


deletion 


(at 37) 


59-62 


elastase 




elastase 




72-80 


deletion 


(at 76) 


elastase 




97-101 


trypsin 


95-102 


elastase 




116 


deletion 


(at 116) 


elastase 




124-133 


addition 


(at 131) 


addition 


(at 131) 


146-151 


elastase 




elastase 




166-179 


trypsin 


164-180 


trypsin 


164-180 


185-187 


trypsin 


184-188 


trypsin 


184-188 


203-206 


elastase 




elastase 




217-224 


chymotrypsin 


216-226 


chymotrypsin 


216-226 



P4 were determined by modeling the substrate on the respective enzyme as 
previously described.^ 

Surface Representations and Electrostatic Potential Maps 

The solvent exclusion surfaces of the substrates and enzymes were calcu- 
lated by the method of Connolly^^'^^ using a program obtained from him. 

The electrostatic potential maps were calculated as follows: a point unit 
positive charge was placed on the solvent-accessible surface and the elec- 
trostatic potential calculated from Coulomb's Law as 

i n 

where is the charge on the ith atom and is the distance from the point 
positive charge to this atom. The charges on the atoms were taken from 
Kollman and his co-workers. A unit dielectric constant is used here because 
we are interested in the interactions between enzyme and substrate at the com- 
plex interface from which water should be excluded. The major conclusions 
of this work are not changed if a dielectric constant = is used as reported 
by many workers. ^^'^^ 

RESULTS AND DISCUSSION 

The Structures of the Enzyme-Substrate Complexes 

The details of the proposed enzyme-substrate complex for factor Xa and 
prothrombin have previously been described^ (Fig. 2). The cleaved peptide 
bond lies between residues Pi and Pi' in prothrombin (see Table 3). Residue 
Pi is an Arg and forms a salt bridge with Asp 189 of factor Xa. This interac- 
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TABLE 3. Sequence of Cleaved Peptide in Substrate 



Substrate 


P4'' 
12* 


Pa 
13 


P2 
14 


P, 
15 


P>' 
16 


P2' 
17 


P3' 
18 


Prothrombin 
Factor X 
Trypsinogen 


He 

Gin 

Asp 


Glu 
Val 
Asp 


Gly 
Val 
Asp 


Arg 
Arg 
Lys 


He 
He 
He 


Val 
Val 
Val 


Glu 
Gly 
Gly 



" Nomenclature as in ref. 17. 

* This residue numbering corresponds to that for chymotrypsinogen which is often used 
tor comparmg serine protease sequences. These numbers are used to label the residues of the 
substrate peptides in Figs. 2-9. 



tion IS typical of many serine proteases and is the molecular basis of the pri- 
mary specificity shown by these enzymes, including trypsin, to cleave only 
after Lys or Arg residues.^' The main chain of residues and P3 forms an 
antiparallel p-sheet with residues 216-218 of the factor Xa enzyme This fea- 
ture appears to be common to all the inhibitor structures examined experimen- 
tally and IS assumed in all the complexes modeled in this study. Residue 
is a Gly and thus has no side chain; any larger side chain at this point would 
collide with the phenolic hydroxyl of Tyr 99 on factor Xa. The glutamate 
side chains at P3 and P/ form specific salt bridges with Arg 143 and Lys 
62 on factor Xa, respectively. Finally, He P^, He P/. and Val P^' all lie in 
hydrophobic regions of the molecule. Thus, to sum up, in addition to the 
Arg-Asp salt bridge conferring the primary specificity, there are two other 
salt bridges, a stereogeometric requirement for a Gly, and numerous 
hydrophobic interactions. 

The complex between blood clotting factor X as a substrate and its ac- 
tivating enzyme factor IXa was also examined in detail (Fig 3) The pri- 
mary specificity interaction is the same as above: Arg P, of factor X with 
Asp 189 of factor IXa. The side chains of Val P, and Val P3 appear to make 
tew close interactions, nor do they seem to lie in a particularly hydrophobic 
environment. On the other hand, Gin P4 has the possibility of making several 
hydrogen bonds with the hydroxyl side chains of Thr 172 and Ser 175 On 
the other side of the cleaved bond, the relatively invariant He P,' and Val 
P2 he m conserved hydrophobic pockets on factor IXa. Residue P3' lies im- 
mediately adjacent to the VR at positions 36-38, which cannot be modeled 
in factor IXa without detailed energetic analysis since no experimentally 
known conformation exists for this sequence (see Methods and Table 2) 
However, the occurrence of Gly P3' in factor X avoids the need to model 
this VR accurately because it has no side chain to fit to the enzyme as there 
was in the prothrombin-factor Xa complex (see above and ref. 7 for a detailed 
discussion). Therefore, the major interactions between factor X and factor 
IXa m addition to the Arg P,-Asp 189 salt bridge, are two or three hydrogen 
bonds and several hydrophobic contacts. 

The contrast between the two enzyme-substrate complexes described above 
IS quite stnking. Whereas the occurrence of two additional specific salt bridges 
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beyond the primary specificity interaction seems likely in factor Xa- 
prothrombin, no such additional charge-charge interactions occur in the 
factor IXa-factor X complex. In fact, the side chains of Val Pa, Val P3', and 
of course, Gly P'3 in factor X seem to contribute very little to either speci- 
ficity or binding to the factor IXa enzyme. 

The very different nature of the interactions found in the factor Xa- 
prothrombin and factor IXa-factor X complexes encouraged us to examine 
another enzyme-substrate interaction in order to see how different these sub- 
strates are from each other. Therefore, the trypsinogen activation sequence 
(Table 3) was examined in the appropriate conformation for binding to an 
activating serine protease, in this case taken to be trypsin (see Fig. 4 and 
Table 1). Once again, the primary specificity salt bridge appears, this time 
between Lys Pi and Asp 189 of trypsin. The three residues prior to the Lys 
are aspartates. No countercharge occurs on trypsin for either Asp^ P2 or 
Asp P3. The £-amino group of Lys 224 on trypsin is approximately 8 A from 
Asp P4, but can be rotated to lie in close association with that carboxylate. 
On the other side of Lys Pi, He P/ and Val P2' lie in their usual, conserved 
hydrophobic environments. Gly at P3', of course, has no side chain to in- 
teract with the enzyme. 

Overall, the trypsinogen peptide seems to interact even less with the trypsin 
active site region than does factor X with factor IXa. In that latter case, H- 
bonds were formed to Gin P4 and valines P2 and P3 did provide some 
hydrophobic interactions beyond the He P/, Val P2' contacts. In the trypsin- 
trypsinogen complex, there is almost no interaction of side chains other than 
that of He P/ and Val P2'. 

This result accords well with the fact that trypsin has lower secondary 
specificity than the other enzymes studied here. The absence of specific, lim- 
iting interactions between the enzyme and substrate side chains at positions 
P2 through P4 and P3' should allow a wide variety of peptides to bind to the 
active site. Since trypsin is presumably intended to cleave denatured and par- 
tially fragmented proteins, any charged or polar side chains on the substrate 
which would not interact with the enzyme (such as Asp P2 through P4) could 
be satisfied by the buffer and the solvent. 

Comparison of Substrate Properties 

It is clear, from examining the three peptide sequences in Table 3 and 
Figs. 2-4, that each substrate will make quite different hydrophilic and 
hydrophobic interactions with its respective cleaving enzyme. Consequently, 



> 

FIGURE 5. Electrostatic potential surface maps for the three substrate peptides of Table 1 
shown in stereo. The electrostatic potential is represented as follows: <- 15 kcal/mole, red; 
- 15 to -5 kcal/mole, orange; -5 to 5 kcal/mole, green; 5 to 15 kcal/mole, light blue; >15 
kcal/mole, blue. Top is prothrombin activation peptide, middle is factor X activation peptide, 
and bottom is trypsinogen activation peptide. The statistics for these maps are presented in Table 4. 
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TABLE 4. Distribution of Electrostatic Potential 
on the Surface Map of the Substrate 


Percent of surface 








area between 


Prothrombin 


Factor X 


Trypsinogen 


< - 15 kcal/mole 


76-.6% 


0.0% 


95.7% 


- 15 to -5 kcal/mole 


8.1 


0.0 


1.4 


- 5 to 5 kcal/mole 


5.1 


0.0 


1.8 


5 to 15 kcal/mole 


3.1 


3.5 


1.1 


> 15 kcal/mole 


7.1 


96.5 


0.0 


Average potential 








(kcal/mole) 


-39.5 


41.3 


-88.0 


Total charge on 








substrate peptide. 


-1 


+ 1 


-2 



it would be useful to compare the different complexes systematically to better 
understand the nature of their specificity. 

Langridge and co-workers^^ have introduced a surface map representa- 
tion of the electrostatic potential of these molecules as a way of depicting 
the charge-charge interactions between molecules in a complex. When this 
method is applied to the substrate peptides (Fig. 5) the differences between 
the various substrates are highlighted. Table 4 gives a breakdown of the per- 





FIGURE 6. A slice through the electrostatic potential surface maps for prothrombin bound 
to the factor Xa enzyme. The positive surface (blue) on the right (labeled R 15) is part of the 
substrate surface at Arg Pi . It is enveloped by a negative (red) surface which is due to the proximity 
to Asp 189 of factor Xa. On the left, the negative surface is that of Glu P3' (labeled E 18). 
Adjacent to it is the positive surface of the enzyme due to Lys 62. Nevertheless, it is very difficult to 
distinguish which surface comes from the substrate and which from the enzyme, without using 
different colors for the two surfaces (rather than colors for the electrostatic potential). 
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cent surface area on each substrate which is positive and negative in order 
to provide a quantitative indication of the differences between these peptides. 

When the electrostatic potential maps are examined in detail, it emerges 
that the charged side chains, Arg, Lys, Glu, and Asp, dominate the poten- 
tial function. Therefore, the map for the prothrombin peptide is positive 
around the Arg P, side chain but negative everywhere else because of the 
two carboxylates of Glu P3 and P3' which give the peptide a net charge of 
- 1 (Table 4). Similarly, the map for the factor X peptide is almost entirely 
positive because the only charged species, the positive Arg at P,, is "felt" 
everywhere on the surface of the peptide. In the same way, the trypsinogen 
peptide with a net charge of -2 due to the three negative Asp residues at 
P2 to P4 overwhelms the effect of the positive Lys at P, except in the im- 
mediate environment of the £-amino group. 



Variation of Substrate-Peptide Electrostatic Potential by Environment 

In order to better understand the implications of the electrostatic poten- 
tial functions, as represented by these color-coded surface maps, to the spec- 
ificity of the substrate-enzyme interaction, the surface maps of several of 
the enzymes were examined in the complex with the respective substrate. For 
this study, the surface maps of the complex of prothrombin with its activating 
enzyme factor Xa (Fig. 2) are shown in Fig. 6. Unfortunately, the superpo- 
sition of the two color-coded surfaces makes them difficult to distinguish; 
consequently, only a thin section can be shown here." It can be seen that, 
for example, the electrostatic potential on the enzyme surface is opposite to 
that on the substrate in the region of the salt bridge interactions described 
above, as shown for Arg P,-Asp 189 and for Glu Pa'-Lys 62. 

Does this mean that the enzyme provides the correct countercharges to 
balance the charged species on the substrate peptide and yield potential values 
close to 0? To test this, the electrostatic potential map at the surface of the 
prothrombin substrate peptide was recomputed, but with the three counter- 
charge residues. Asp 189, Arg 143, and Lys 62 included to give a total charge 
of zero for the substrate + countercharges system. The resulting map is shown 
in Fig. 7B and summarized in Table 5. The charge distribution is now strik- 
ingly different than for the substrate peptide by itself (Fig. 7A). The first 
point to stress is that very little of the surface area, 13%, is close to 0 poten- 
tial. Thus, the countercharges are not neutralizing the charge on the sub- 
strate but are changing the electrostatic potential distribution considerably. 
The surface about P,' to P3' changes from negative to positive because of 
its proximity to Lys 62, Arg P^, and Arg 143. Similarly, the tip of the sur- 
face for Arg P, becomes negative because of the close approach of Asp 189. 
Some of the molecular surface does reflect the nature of the atoms that are 
immediately nearby, rather than longer range charge effects, as for example, 
the side chain of He P4 is close to zero whereas its carbonyl oxygen is negative. 

Having seen the effect of including just the countercharges, how is the 
electrostatic potential of substrate peptide affected by being placed in the 
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TABLE 5. Distribution of Electrostatic Potential 



for the Prothrombin Substrate Peptide 


Percf*nf nf siirfarf* 




rruinroiTiDin + 


rroinromDin + 


Prothrombin + 


area between 


Prothrombin 


3 countercharges 


factor Xa 


factor Xa + Ca* 


< - 15 kcal/mole 


76.6% 


24.4 <7o 


99.9% 


70.1<7o 


- 15 to -5 kcal/mole 


8.1 


11.1 


0.1 


13.3 


- 5 to 5 kcal/mole 


5.1 


13.1 


0.0 


8.2 


5 to 15 kcal/mole 


3.1 


16.5 


0.0 


6.5 


> 15 kcal/mole 


7.1 


34.9 


0.0 


1.9 


Average potential 










(kcal/mole) 


-39.5 


2.6 


-58.4 


-27.3 


Net charge on system 


-1 


0 


-1 


+ 1 



complete environment of the factor Xa enzyme model structure? The net 
charge of the factor Xa molecule is 0. This gives a total charge of - 1 for 
the enzyme-substrate complex. The result of this calculation was a complete 
surprise; Fig. 7C and Table 5 present the resuhs. The entire map is negative 
with an average electrostatic potential of - 58 kcal/mole. Thus, the peptide 
in the active site of the enzyme must lie in a highly negative environment 

< — 

FIGURE 7. Electrostatic potential surface maps for the substrate prothrombin activation peptide 
(see Table 5). (A) Potential map for prothrombin peptide by itself with a net charge of - 1. 
(B) Map for the prothrombin peptide with the three salt bridge residues: Lys 62, Arg 143, and 
Asp 189. The net charge on this system is 0. (C) Potential surface map for the prothrombin 
peptide.on the factor Xa enzyme. Note that the whole surface is now highly negative with mean 
potential value of -58 kcal/mole. This is due to the distribution of charged residues on the 
enzyme (see Table 6 and text for discussion). (D) Ca** placed 14 A from the substrate near 
residues Glu 36 and Glu 37 of factor Xa (see Fig. 8), Note that the positive surface corresponds 
exactly to the carboxylate side chain of Glu P3' (labeled E 18) which is itself negative. 

FIGURE 9. Electrostatic potential surface maps for the substrate factor X activation peptide 
(see Table 7). (A) The map for factor X peptide which has a net charge of + 1, The map is 
dominated by this positive charge of Arg P,. (B) Potential map for the peptide together with 
Asp 189, the countercharge for Arg Pi (labeled R 15). The net charge is now 0. A discussion 
of this distribution appears in the text. (C) Map for the factor X peptide together with the 
model for its activating enzyme, factor IXa. In the factor IXa enzyme, the active site and speci- 
ficity sites are highly positive, in contrast to factor Xa (FiG.7C), where the environment was 
highly negative. The potential at the surface of residue Val P2' (labeled V 17) is negative, due 
to a concentration of negatively charged residues in the region of the enzyme adjacent to this 
position (see text). 

FIGURE 10. Electrostatic potential surface maps for the substrate trypsinogen activation peptide 
(see Table 8). (A) The map for the trypsinogen peptide alone with a charge of -2 due to the 
Asp residues. The surface is overwhelmingly negative, (B) The trypsinogen peptide in the environ- 
ment of trypsin gives this potential map. The A^-terminal portion remains negative since there 
are no countercharges for Asp P2 and Asp P3 (labeled D 14 and D 13, respectively) and they 
point out into solvent. There is a countercharge for Asp P4 (labeled D 12) which is Lys 224, 
and thus some positive surface appears at this residue. The C-terminal portion is apparently 
dominated by the large positive net charge on the trypsin molecule of + 10. Note that the transition 
from negative to positive occurs right at the point of the cleaved peptide bond. 
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even though, as noted above, the net charge on the enzyme is 0. In order 
to determine why this is the case, the net charge on the enzyme was calcu- 
lated in shells around the substrate, as shown in Table 6. As can be seen, 
factor Xa has an excess of negative charges within 7 A of the substrate. There 
are shells with an excess of positive charge also, but these lie much farther 
away, from 12 to 20 A, and thus by Coulomb's law, have a much smaller 
effect on the substrate. It is particularly striking that the active-site environ- 
ment of factor Xa is so negative, when the net charge on the prothrombin 
peptide is also negative. 

Is this highly negative environment of the active site a realistic view of 
this enzyme-substrate interaction in the factor Xa-prothrombin complex? 
Clearly, these charge effects could be modified and modulated by the binding 
of ions from solvent and buffer onto the protein. In particular, there is a 
loop in factor Xa that contains three Glu residues very close together, at po- 
sitions 36, 37 and 39 (Fig. 8). Two, and perhaps all three, of these residues 
may form a Ca** ion binding site. In order to test the effect of the binding 
of such an ion, a Ca** was bound to the enzyme at this site and given a charge 
of + 2.0. It is important to point out that the closest this ion approaches 
to the substrate peptide is 14.1 A at Glu P3' (see Fig. 8). 

The electrostatic potential map of the prothrombin activation peptide was 
recalculated in the environment of the factor Xa enzyme model with the one 
Ca** ion bound, as described. The resulting map is shown in Fig. 7D and 
is summarized in Table 5. There is a significant influence of this single Ca** 
ion on the electrostatic potential in the environment of Glu P3' and its sur- 
rounding, which has changed from significantly negative to positive. Thus, 
it is clear that a proper treatment of the electrostatic environment of the sub- 
strate must take into account solvent and buffer ion interaction with the en- 
zyme even at sites somewhat removed from the immediate region of the ac- 
tive and specificity sites of the enzyme. 

Similar calculations were performed on the electrostatic potential map 
of the factor X activation peptide substrate bound to its activating enzyme 
factor IXa (see Fig. 9 and Table 7). As previously discussed, the substrate 
peptide appears to be entirely positive with a mean electrostatic potential value 
of 41 kcal/mole (Fig. 9A) due to the Arg at the Pi site. When the counter- 
charge for this Arg, Asp 189 from the protein, is included in the calculation 
of the electrostatic potential, the result is a significant difference in the poten- 
tial function (Fig. 9B and Table 7) which now has an average potential value 
of 8.4 kcal/mole. Very little of the surface is actually negative, mostly in 
the immediate region of Asp 189. There is a positive band that is largely due 
to the influence of the Arg P, guanidinium group. Moving further away from 
Arg P„ there is a belt of approximately neutral or slightly positive or nega- 
tive potential that is largely the result of the carbonyl oxygens of residues 
P2 and P4 and the Gin P4 side chain. Finally the edge of the molecule which 
points mostly away from the activating enzyme (not easily seen in this view) 
is again positive. This appears to be due to a series of main chain amino 
hydrogens of residues P2, Pi and P/ all pointing in a similar direction (see 
Fig. 3). Once again the addition of a countercharge, which gives a formally 
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TABLE 7. Distribution of Electrostatic Potential 
for the Factor X Substrate Peptide 



Percent surface 


racior A 


Factor X 
+ 1 countercharge 


Factor X 
+ factor IXa 


< - 15 kcal/mole 
- 15 to -5 kcal/mole 
-5 to 5 kcal/mole 
5 to 15 kcal/mole 
> 15 kcal/mole 


0.0% 
0.0 
0.0 
3.5 
96.5 


4.3% 

7.7 
24.7 
36.1 
27.3 


0.0% 
1.9 
1.9 
3.0 
93.2 


Average potential 
(kcal/mole) 


41.3 


8.4 


64.9 


Net charge on system 


1 


0 


2 



net uncharged system, results in a distribution of very significant positive 
and negative potentials rather than a close to zero value. 

Next, the electrostatic potential map was recalculated for the factor X 
peptide in the environment of the factor IXa model structure (Fig. 9C and 
Table 7). The resulting map is different from both the isolated peptide and 
the peptide-countercharge distribution. The potential map is almost entirely 
positive. This is not the result of the overall + 1 charge on the factor IXa 
molecule (versus - 1 for factor Xa, see Table 6) because the addition of + 2 
to factor Xa (Fig. 7D) had only a limited effect on the negative environment 
even though the total net charge of factor Xa + Ca** is + 1. Rather, it is 
once again the result of the charge distribution on factor IXa (Table 6). 
Whereas the immediate region about the substrate from 0-7 A tends to be 
negative, a larger net positive charge from 7-15 A seems to overwhelm the 
negative group and gives the active site a positive character. 

Interestingly, despite the overall positive environment, the surface about 
the side chain and carbonyl oxygen of Val P^' is quite negative. This is due 
to a concentration of negative charges and a paucity of positive charges in 
this region of the factor IXa enzyme, in particular, Glu 36 and Glu 74 but 
also Glu residues at positions 70, 78, and 80. 



TABLE 8. Distribution of Electrostatic Potential 
for the Trypsinogen Substrate Peptide 



Percent surface 
area between 


Trypsinogen 


Trypsinogen 
+ trypsin 


< - 15 kcal/mole 
- 15 to -5 kcal/mole 
-5 to 5 kcal/mole 
5 to 15 kcal/mole 
> 15 kcal/mole 


95.7% 
1.4 
1.8 
1.1 
0.0 


30.1% 

9.4 

5.9 

3.4 
51.2 


Average potential 
(kcal/mole) 


-88.0 


21.3 


Net charge on system 


-2 


+ 8 
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The average potential of 65 kcal/mole for factor X on the factor IXa 
enzyme indicates that the factor X substrate peptide, which itself has a net 
positive charge, resides in a highly positive environment in factor IXa. This 
is the exact parallel of the situation in the enzyme factor Xa, which was a 
negative substrate peptide in a highly negative environment (see Fig. 7C). 
Thus, in each case, the enzyme active site carries a potential that has the same 
sign as the potential on the substrate peptide. One might speculate that this 
may represent a mechanism whereby the enzyme prevents too strong binding 
of the substrate -which would result in diminished turnover. 

Similarly dramatic differences can be observed in the electrostatic poten- 
tial function of the trypsinogen activation peptide on the enzyme trypsin (Fig. 
10 and Table 8). In the isolated substrate peptide, as noted above, the three 
Asp residues at Pa through P4 dominate the potential function, virtually over- 
whelming the positive charge at Lys Because countercharges for Asp P^, 
Asp P3 and Asp P4 (except perhaps Lys 224 for the latter, see Fig. 4), do 
not occur on trypsin, the potential map for the substrate peptide together 
with enzyme countercharges cannot be computed. However, the addition of 
the enzyme which has a net -1- 10 charge has a remarkable effect on the poten- 
tial function. The A^-terminal half of the substrate peptide remains negative, 
due presumably to the aspartate residues. Only Asp P4 has a possible coun- 
tercharge on the enzyme in residue Lys 224, hence the positive and zero re- 
gion at the surface of residue P4. The C-terminal half of the molecule does 
apparently "feel" the highly positive nature of the molecule. The complete 
substrate is not dominated by the -f- 10 charge of trypsin because most of 
the positive side chains of trypsin are distant from the substrate (see Table 
6). It is especially clear that both with respect to the Asp residues at posi- 
tions P2 through P4 and to the many excess positive residues on trypsin, all 
of which are exposed to solvent, the binding of ions from the buffer and 
solvent should play a crucial role in the proper evaluation of the electrostatics 
of substrate-enzyme interaction. 

Because trypsin is a digestive enzyme, it tends to be less specific than the 
other serine proteases studied here. Therefore, it probably operates mostly 
on denatured proteins or highly accessible peptide portions of macro- 
molecules. Thus, it is likely that a substrate peptide such as in trypsinogen 
will bind to trypsin with Asp residues P^ and P3 and perhaps even P4 (de- 
spite the Lys 224 counterion) interacting ultimately with solvent and buffer 
ions. Hence, the electrostatic effect of these aspartates would be consider- 
ably diminished. It will be interesting to examine several other possible sub- 
strates of trypsin in order to see how they are affected by the trypsin envi- 
ronment and compare these results with those obtained for the more specific 
blood clotting factors described above. 
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